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ION HEATING IN A PLASMA FOCUS 

By Frank Hohl and S. Peter Gary* 

Langley Research Center 

SUMMARY 

Ion acceleration and heating in a plasma focus are investigated by the numerical 
integration of the three-dimensional equations of motion. The electric and magnetic fields 
given were derived from experimental data. The results obtained show that during the 
collapse phase of focus formation, ions are efficiently heated to temperatures of several 
keV. During the phase of rapid current reduction, ions are accelerated to large velocities 
in the axial direction. The results obtained with the model are in general agreement with 
experimental results. 

INTRODUCTION 

Plasma focus devices (ref. 1) are presently under study in a number of laboratories. 
When operated with deuterium at pressures of a few hundred pascals (a few torr), the 
plasma focus emits neutron bursts of 10^ to 10l3 neutrons with a neutron production time 
around 100 ns. Average energy shifts of up to 500 keV, corresponding to center-of-mass 
velocities up to 2 Mm/ s are obtained by analyzing the neutron energies in the axial direc- 
tion (refs. 2 and 3). In general, the experimental data obtained by the various investiga- 
tors (refs. 2 to 8) do not fit either of the simple models commonly proposed, namely, the 
beam-target or moving-boiler models. The experimental observations which a successful 
model of the plasma focus formation must explain are ion heating to keV temperatures 
(refs. 3 and 9) and axial ion velocities corresponding to energies of at least 300 keV 
(refs. 2, 10, and 11). 

The two-dimensional fluid model of Potter (ref. 12) describes well the initial phase 
of focus formation and agrees with a number of the observed properties. Such a model, 
however, does not include the detailed process of ion and electron heating apparently 
occurring during the final phase of focus formation. For example, the two-dimensional 
magnetohydrodynamic simulations are unable to account for the experimentally observed 
anisotropy in the neutron emission. 

Bernstein (ref. 13) has presented a model that more closely represents the experi- 
mental data. The present work extends his results but differs in a number of important 
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aspects. First, in the present work two variations in the electric field are considered: 
one resulting from the radial velocity of the collapsing plasma and the other resulting 
from the rapid current decrease during the final focus formation. Bernstein considered 
only the electric field induced during the collapsing phase and thus, neglected the important 
phase of rapid current reduction. Experimentally, this current- reduction phase has been 
proven to be a characteristic of focus formation. 

Also, in his model Bernstein assumed an on-axis electric field of zero; this assump- 
tion implies a strong electric field away from the axis, where the magnetic field is strong. 
This leads to heating, but also prevents ions from acquiring large axial velocities. In the 
present work, the axial electric field is zero outside the plasma (at the outer electrode), 
in agreement with Lenz' law. Thus, the strong electric field in the region of low magnetic 
field near the axis allows strong ion acceleration. Since experimental results indicate a 
correlation between the rapid current drop and ion acceleration, this phase of focus forma- 
tion is also investigated in the present work. The present work improves Bernstein's 
work in several other respects in that it calculates the ion motion in three-dimensional 
space rather than in two, treats initial distribution with Maxwellian velocity distributions, 
and considers much larger numbers of ions. 

It should be emphasized that the present calculation, like Bernstein's, is not self- 
consistent. The electric and magnetic fields are assumed as given, and the response of 
the ions to the fields is calculated. Assumii^ the fields as given is possible only when 
the ion response does not significantly affect the fields. Since the electrons are the pri- 
mary current carriers, this assumption is satisfied for the present work. The present 
method has the advantage of allowing a fully three-dimensional treatment that includes 
the effects of strongly inhomogeneous fields. Some initial results of the present work 
were published elsewhere (ref. 14). 

Figure 1 displays a cross section through a plasma focus apparatus, which consists 
of two coaxial cylindrical electrodes. The outer electrode has an inner diameter of 10 cm 
and the inner electrode has an outer diameter of 5 cm. Figure 2 shows a streak photo- 
graph of the plasma focus formation taken with an image converter camera. The photo- 
graph displays a region 0.5 cm above the center electrode viewed through a 0.2-cm-wide 
slit perpendicular to the axis. The time interval from -50 ns to 0, as shown in figure 2, 
when the plasma column is in the final collapse phase, is considered the compression 
phase in the present paper. The time interval from 0 to about 50 ns is the period when 
maximum compression has occurred, and during this time, a rapid current reduction takes 
place. This rapid current reduction is illustrated in figure 3, which shows typical oscillo- 
scope traces of the current and induced electric field in the plasma focus apparatus. 
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Note that because of the relatively slow detector response, the current reduction appears 
to take place during a time of about 200 ns. However, the electric field induced by the 
changing current shows that the principal change in the current takes place during a time 
of the order of 50 ns. 

SYMBOLS 

B magnetic field 

E electric field 


e 

I 

j 

jo 

h 

m 

r,0,z 

rc 

^0 

^1 

^2 

T 

T* 

t 


magnitude of electronic charge 
total current 
current density 

current density for the compression phase at t = 0 
current density for the rapid current- reduction phase 
mass of deuterium ion 
cylindrical coordinates 

radius of current-carrying region at collapse 
initial radius of plasma column for the compression phase 
radius of plasma column at maximum compression 
radius of outer electrode 


temperature derived from ^ m^V - 
effective energy of ion beam derived from -i m 


time 
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V 


velocity 


Vf, collapse velocity 

x,y,z rectangular coordinates 

|3 inverse time constant 

permeability of free space 

Subscripts: 

max maximum 

x,y,z X-, y-, and z-components 

0 azimuthal component 

Notation: 

( ) root- mean- square value 

unit vector 


THE MODEL 


The equations describing the motion of deuterium ions subject to electric and mag- 
netic fields are 


^ = -^[E(r,t) + V(t) X B(r,t)] 


( 1 ) 


and 


5i“ = v(t) 

dt 


( 2 ) 


where cylindrical symmetry of the fields is assumed. The electric and magnetic fields 
are determined from Maxwell’s equations 


V X E = - — 

at 


( 3 ) 
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and 


V X B = MqT (4) 

where Mq is the permeability of free space, J is the current density, and the displace- 
ment current is neglected. The current density is assumed to be azimuthally symmetric 
along the axis of the plasma focus; that is 

T(r,t) = (5) 

Thus, 

B(r,t) = 0B0(r,t) (6) 

and 

E(r,t) = zEz(r,t) ( 7 ) 

Two phases of the current variation in the plasma focus are investigated. The first 

is that occurring during the compression phase when a circular region of uniform current 
density is compressed. Nonuniform current densities are not expected to change signifi- 
cantly the results obtained (ref. 13). The second phase of current variation is that which 
occurs near peak compression and results in a rapid reduction of the total current. 

Compression Phase 

The final compression phase of plasma focus formation occurs over such a short 
time interval that the total current can be assumed to remain constant (refs. 1 and 8). 

Also, the current distribution is assumed to be uniform over a circular region, the edge 
of which is compressed at a radial velocity V^, of 180 km/s (ref. 15). The radius of the 
circular region is given by 


^c = 


Vet 


( 8 ) 


where ro is the initial radius at t = 0. Thus, the axial current density distribution is 


r 2 

iz(r»t) = io ^ (r ^ rj.(t)) 

iz(r,t) = 0 (r > rc(t)) 


( 9 ) 
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The resulting magnetic field is 

2 r 2 

c 

2 r2 

From Maxwell's equation 




(r ^ r,(t)) 
(r > rc(t)) ^ 


( 10 ) 


9Ez _ 

ar ~ at 

the axial component of the electric field is found to be 


E 


:(r,tl . (r^ - r,2) (r S r,(t))' 

2r„'^ 


E„(r,t) = 0 


(r > rc(t)) 


( 11 ) 


Rapid Current- Reduction Phase 

Near the time of maximum compression, the characteristic current trace for the 
plasma focus shows a reduction in the total current of the order of 30 percent during a 
time period about 50 ns. Thus, the current density variation during this time period is 


iz(r,t) = jjd - ^t) 

j2(r,t) = 0 


(r > r,) 


( 12 ) 


where r, is the radius of the plasma column at maximum compression, and )j = l/ar,^, 
where I is the total current through the focus. This current distribution results in a 
magnetic field given by 


B0(r,t) = (1 - 0t) 

B0(r,t)=-^^^^(l-^t) 


) 

(r2>r>r^) 


(13) 
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The corresponding electric field is 



j 


The radius r 2 > which corresponds to the outer electrode of the plasma focus, is 5 cm for 
the present investigation. The induced electric field outside rg is zero, and this bound- 
ary condition is used in deriving equation (14). Note that the electric field given by that 
equation is time independent, since the current decrease is linear. 


RESULTS 


The first case considered is the collapse phase of the plasma focus. The magnetic 
and induced electric fields are given by equations (10) and (11), respectively. The initial 
radius of the collapsing plasma column r^ is taken to be 1 cm and the collapse speed 
of the column Vc is 180 km/s (ref. 15). Also, the total current is 1 MA, so that 

j = ^ (12 GA/m^j in equations (10) and (11). At the beginning of the simulation, the 

TjTo^ \7T / 

deuterium ions are uniformly distributed over a circular region of radius rQ and have 
a Maxwellian velocity distribution corresponding to a temperature of 10 eV (ref. 16). 
Results obtained with an initial temperature of 100 eV are essentially identical with those 
presented herein. 

Figure 4 shows the variation in the magnetic field at six times during the collapse 
phase. The magnetic field near the edge of the current-carrying column increases to 
large values as the column is compressed. Figure 5 shows the corresponding variation 
in the induced electric field. As can be seen from these two figures, the magnitude of the 
electric field peaks at the axis (r = 0), where the magnetic field is zero. However, since 
the E X B ion drift is radially outward, the condition of a large axial electric, field and 
a low magnetic field cannot be effectively exploited to give large axial ion velocities. 

This is illustrated in figure 6, which displays the evolution of the spatial ion distribution 
and shows that the ions are expelled from the central region by the E x B drift. Never- 
theless, there occurs strong heating for an appreciable portion of the ions, as illustrated 
in figure 7, which shows the evolution of the transverse ion velocities in Vx~Vy space. 

The evolution of the axial velocities as a function of axial position is shown in fig- 
ure 8. Again, large velocities are acquired by a considerable fraction of the ions, and 
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at the same time, the ions perform a cyclotron motion which confines them near their 
original location of z = 0. These results are summarized in figure 9, which displays the 
evolution of various temperature components derived from the rms velocities. In this 
figure Txy represents the transverse temperature, Tz is the axial component, and 
Tjj* corresponds to the energy in directed beam motion in the z-direction. As can be 
seen, a temperature of several keV is obtained as a result of heating due to induced elec- 
tric fields during the collapse phase. Note also that this heating takes place during the 
time that the ions are expelled from the central region. 

As the E X B drift forces the ions from the central region, the mean radius <r> 
of the ion distribution increases from about 7 mm to 8 mm. After 20 ns, the mean radius 
and the temperature show little further increase. When all the ions are initially confined 
to a smaller circular region of radius 1 mm, then the mean radius <r> increases by a 
factor of 10, from about 0.7 mm to 7 mm. This is illustrated in figure 10, which shows 
the evolution of a system identical with that in figure 6, except that now the initial radius 
of the uniform ion distribution is reduced by a factor of 10. The temperature components 
calculated for the system in figure 10 are also larger by a factor of 10 than those shown 
in figure 9. Again, the heating occurs primarily during the initial 20 ns as the ions are 
expelled from the central region. 

The results for the rapid current- reduction phase are presented next. Equa- 
tions (13) and (14) now give the magnetic and induced electric fields. The current reduc- 
tion takes place near the time of maximum compression, when the radius of the current- 
carrying column is about 1 mm. Thus, where rj = 1 mm and the total 

current I is 1 MA. The parameter )3 in equation (12), which is determined in agree- 
ment with experimental data to obtain a 30-percent current reduction in 50 ns is 1 Ms”^. 
The radius Y 2 in equation (14) corresponds to the radius of the outer electrode of the 
plasma focus apparatus and has a value of 5 cm. 

The variation in the magnetic field, as given by equation (13), is shown in figure 11 
at six times during the current- reduction phase. Figure 12 shows the induced electric 
field, which remains constant during the current-reduction phase, and figure 13 shows the 
evolution of the initially uniform ion distribution in x-y space. Since the compression 
phase, with its strong ion heating, precedes the phase of rapid current decrease, a larger 
initial temperature of 200 eV was given to the ions. The initial radius of the ion distri- 
bution in figure 13 is 1 cm. Note that the E x B drift is now radially inward and causes 
an inward compression of the ions. 

In the central region of large electric field and low magnetic field, the ions are 
rapidly accelerated along the z-axis. At the same time, more ions are pushed toward the 
axis by the E x B drift. Figure 14 summarizes some of these results. The transverse 
temperature increases by only about a factor of 4, while Tz shows a much larger 
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increase because new ions continually reach the central region, where they are accel- 
erated along the z-axis, thus producing a large axial velocity spread. 

Since most of the acceleration occurs near the axis, an initial ion distribution iden- 
tical with that shown in figure 13, except that the initial radius was reduced to 1 mm, was 
investigated. The resulting evolution of the spatial ion distribution is shown in figure 15. 
Note that during the time interval of 50 ns shown, all the ions are compressed to a small 
radius. Figure 16 shows a typical ion orbit in both the x-y and the r-z planes. 

Starting at a radius of about 0.6 mm, the ion begins to spiral inward toward the axis. 

After reaching a radius of about 0.2 mm, the ion begins to circle the z-axis and at the 
same time experiences a large acceleration along the z-axis. 

Figure 17 shows the evolution of the Vy-Vjj distribution of the ions. Since all the 
ions are now in the region near the axis, they experience a considerable increase in tem- 
perature, as shown in figure 18, where these results are summarized. First, it can be 
seen that the transverse energy of the ions near the axis reaches several keV. At the 
same time, Tz reaches a value of about 60 keV. The value of Tz shown in figure 18 
should not be interpreted as a true temperature, as it results primarily from the different 
times at which ions reach the central region, where efficient acceleration along the z-axis 
occurs. Note that after all ions have arrived in this region (t « 20 ns), Tz no longer 
increases. The average beam energy of the ions Tz* reaches very large values in this 
simulation because all the ions now acquire large velocities along the z-axis. 

Since the accelerating region, or current column, in the plasma focus extends only 
about 3 cm, the acceleration of ions should terminate when z is greater than that value. 
Figure 19 shows that the fastest ions reach z = 3 cm at about t = 30 ns, at which time 
Tz* is approximately 500 keV (fig. 18). In the actual plasma focus, Tz* should level 
off near that value. 

CONCLUDING REMARKS 

The model presented provides for both ion heating to several keV and ion accelera- 
tion in the axial direction to velocities corresponding to an energy of several hundred keV, 
in agreement with experimental results. The initial collapse phase is the primary ion 
heating phase, whereas the high axial velocity of the ions is acquired during the phase of 
rapid current reduction. The fields were determined from experimental data, and the ions 
moved under the influence of these given fields. Taking the fields as given is justified 
since the much less massive electrons are the primary current carriers and thus deter- 
mine the structure of the magnetic and induced electric fields.' A more detailed under- 
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standing of focus formation, such as the cause of the increased resistivity which produces 
the rapid current drop, will require more involved simulations and analyses. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 12, 1974. 
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Figure 1.- Cross section through coaxial plasma focus apparatus. 
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Figure 2.- Streak photograph of focus formation. (The photograph 
displays a region 0.5 cm above the inner electrode, viewed 
through a 0.2-cm-wide slit perpendicular to the axis.) 




Figure 3.- Traces of typical oscillograms of current and induced electric field. 
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Figure 4.- Radial variation in magnetic field at six times 
during collapse phase of focus formation. 


t = 20 ns t = to ns 



Figure 5.- Radial variation in induced axial electric field at six times during 
collapse phase of focus formation. (Note that the electric field peaks on 
the axis and reaches values in excess of 10 MV/ m.) 
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Figure 6.- Evolution of ion distribution during collapse phase of focus formation. 
(Note that ions are expelled from the central region by the E x B drift.) 
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Figure 7.- Evolution of transverse ion velocity distribution during collapse phase. 
(A fraction of the ions are heated to keV temperatures.) 
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Figure 8.- Evolution of ion distribution in Z-V 2 space. (Note that large axial 
velocities are acquired by a considerable fraction of' the ions. At the same 
time the ions are confined near z = 0 because of their cyclotron motion.) 
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Figure 11.- Radial variation in magnetic field at six times 
during current- reduction phase of focus formation. 



Figure 12.- Radial variation in axial electric field during current- reduction phase. 
(Note the strong electric field of several megavolts per meter on the axis.) 
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Figure 15.- Evolution of central ion distribution during current- reduction phase. (Note that 
the initial radius is 1 mm, or one-tenth the radius of the distribution in fig. 13.) 
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jectory in x-y space. 


t = 10 ns 

1 


(b) Trajectory in r-z space. 

Figure 16.- Typical ion trajectory during current- reduction phase. (The ion spirals 
toward the axis. As it nears the axis it begins to circle the axis and at the 
same time is accelerated in the z-direction.) 
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Figure 17,- Evolution of ion velocity distribution during current- reduction phase. 
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